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ABSTRACT 

We present a technique that appHes spectral synthesis to medium resohition spectroscopy (MRS, 
R ^ 6000) in the red (6300 A < A < 9100 A) to measure [Fe/H] and [a/Fe] of individual red giant stars 
over a wide metallicity range. We apply our technique to 264 red giant stars in seven Galactic globular 
clusters and demonstrate that it reproduces the metallicities and a enhancements derived from high 
resolution spectroscopy (HRS). The MRS technique excludes the three Ca II triplet lines and instead 
relies on a plethora of weaker lines. Unlike empirical metallicity estimators, such as the equivalent 
width of the Ca II triplet, the synthetic method presented here is applicable over an arbitrarily wide 
metallicity range and is independent of assumptions about the a enhancement. Estimates of cluster 
mean [Fe/H] from different HRS studies show typical scatter of ~ 0.1 dex but can be larger than 
0.2 dex for metal-rich clusters. The scatter in HRS abundance estimates among individual stars in a 
given cluster is also comparable to 0.1 dex. By comparison, the scatter among MRS [Fc/H] estimates 
of individual stars in a given cluster is ~ 0.1 dex for most clusters but 0.17 dex for the most metal-rich 
cluster, M71 (([Fe/H]) = —0.8). A star-by-star comparison of HRS vs. MRS [a/Fe] estimates indicates 
that the precision in [a/Fe]j^j^g is 0.05 dex. The errors in [Fe/H]]yjj^g and [a/Fe]]y[j^g increase beyond 

0.25 dex only below signal-to- noise ratios of 20 A , which is typical for existing MRS of the red giant 
stars in Leo I, one of the most distant Milky Way satellites (250 kpc). 

Subject headings: stars: abundances — globular clusters: individual (M79, NGC 2419, M13, M71, 
NGC 7006, M15, NGC 7492) 



1. INTRODUCTION 

In the popularly accepted paradigm of hierarchical 
structure formation (Searle fc Z inn 1978; IWhite fc ReesI 
Il978[ ) , large galaxies grow by assembling smaller compo- 
nents. Simulations motivated by ACDM cosmology have 
successfully reprod uced the detailed properties of ob- 
served galaxies (e.g.. lRobertson et al.ll2005tlDutton et al.l 
[2005). Such simulations have advanced enough to track 
an important prediction of ACDM cosmology: the chemi- 
cal properties of stars in different dynamical components 
in a single dark matter halo, including the cold dwarf 
satellites, the dissol ving tidal st reams, and the hot stellar 
halo. For example, iFont et all (2006) make predictions 
of metallicity and a enhancements of a Milky Way-like 
halo. They verify their metallici ty predictions wi th pho- 
tometric measurements of M31 (jFont et al.ll2008[ ). 

Several techniques can be used to measure stellar abun- 
dances. The most trusted abundance measurement tool 
is high resolution spectroscopy (HRS). A high resolu- 
tion stellar spectrum contains information about the 
temperature, surface gravity, and abundances of indi- 
vidual elements in the stellar atmosphere. However, 
HRS spreads the light of a star over finely spaced res- 
olution elements and therefore requires long exposures 
of bright stars. With the current state-of-the-art tele- 
scopes, typical HRS targets are brighter than V ^ 15, 
although some authors have observed targets as faint as 



y - 18 fe.g.. iCohen fc Melended [2005ah . HRS targets 
have ranged from moderate to low luminosity stars in 
the solar neighborhood to high luminosity stars through- 
out t he Milky Way (MW) a nd even in its dwarf satellites 
(e.g.. IShetrone et al.l[2003[ ). HRS within the MW has 
allowed measurements of [Fe/H], [a/Fe], and individual 
element abundances of in dividual stars in d ifferent dy- 
namical components (e.g., I Venn et al.ll2004h . However, 
it is currently not feasible to obtain large HRS samples 
of individual stars beyond the MW and its satellites. 

In order to reach fainter and more distant stellar sys- 
tems, photometry of resolved stellar populations is a 
commonly used abundance measurement tool. Most 
studies based on color-magnitude diagrams (CMDs) tar- 
get the most luminous components of a population — the 
red giant branch (RGB) and horizontal branch (HB) — 
because they are visible over large distances. CMDs 
that do not reach the main sequence turn-off (MSTO) 
are susceptible to the age-metallicity degeneracy and 
the H B second parameter pro blem (Sandag c fc W ildej 
19671). As a rare exception, iBrown et al.l (|2003t ) and 



Brown et al.l (|2006( l have used extremely deep imaging 



of three fields in M31 to constrain not only metallicity 
distributions but also ages of stellar populations. How- 
ever, such deep photometry of M31 is extremely expen- 
sive and infeasible for more than a few fields. Even with 
the MSTO, photometry provides only one abundance di- 
mension, overall metallicity. In fact, photometric abun- 
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dance estimates based on the RGB must assume [a/Fe] . 

Medium resolution spectroscopy (MRS) avoids most 
of the assumptions involved in photometric metallicity 
estimates. Furthermore, it can provide multiple abun- 
dance dimensions, such as [Fe/H] and [a/Fc]. State- 
of-the-art MRS instruments include lris and deimos 
on the Keck telescopes, fors and VIMOS on the Very 
Large Telescope, and IMACS on Magellan. Some of these 
instruments have been used to obtain MRS of RGB 
stars in the bulge, disk, halo, and dwarf satellites of 
M31 (iGuhathak urta erall l200 6l: IChapman et"all l2006t 
iKalirai et aLl l2006: Koch ct al. 2007t). 

The traditional MRS abunda n ce te chnique is spec- 
trophotometric indices. iPrestonI (jl961h established the 
first calibration between the equivalent width (EW) of 
a single metal line (Ca II K A3933) and stellar metal- 
licity. The Ca II AAA8498, 8542, 8662 infrared triplet 
has been a more popular metal l icity indicator in the 
last decade. lArmandroff fc Ziiml (|1988f ) first quantified 
the correlation between the EW of the Ca II triplet in 
integrated light spectra of MW GCs and their previ- 
ously measured [Fe/H] . A plethora of astronomers have 
since develop ed their own Ca II tri plet metallicity cali- 
bratio ns, but [Olszewski et all ()199lD and lRutledge et all 
()1997( ) are the most highly cited. Despite the success of 
the Ca II triplet as a metallicity indicator over a wide 
range (—2.2 < [Fe/H] < —0.6), its use in determining 
[Fe/H] necessit ates assuming the ratio [Ca/Fe]. In fact, 
iRutledge et al.l warn, "caution — perhaps considerable — 
may be advisable when using [Ca II triplet reduced 
width] as a surrogate for metallicity, especially for sys- 
tems where ranges in age and metallicity are likely." For 
an analysis of the different effects of [Fe/H] and [ Ca/H] 
on the Ca II triplet EW, see lBattaglia et all (|2008[ ). who 
also claim errors of 0.1 dex in [Fe/H] at signal-to- noise 

ratios S/N > 10 A . Finally, any empirical calibration 
is restricted to the metallicity range of the calibrators. 
The Ca II triplet met allicities may not be accurate be- 
low [Fe/H] < -2.2 (K och et al.| [2008). 

Spectral modeling, essentially the HRS abundance 
technique, circumvents dependence on the properties of 
calibrators. The procedure is to synthesize a spectrum of 
a stellar atmosphere and compare to an observed spec- 
trum. The atmospheric parameters and abundances can 
be adjusted until the best fit is achieved. Most HRS stud- 
ies model individual lines and compare synthetic and ob- 
served EWs. At medium and low resolution, most metal 
absorption lines are weak, blended, or both. However, a 
few thousand A of spectral coverage includes hundreds of 
metal absorption lines. Pixel-to-pixcl spectral fits lever- 
age the statistical power of many lines to find the optimal 
atmosphere and abundance. 

Thi s met hod is not new t o MR S. For example. ISuntzef^ 
(jl981h and ICarbon et all ()1982f ) modeled blue spectra 
at low resolution to determine C and N abundances 
in globular cluster (GC) RGB stars. Synthetic pixel- 
to-pixel matching is one metallicity estimator in the 
Sloan Extension for Galactic Understanding a nd Explo- 
ration (S EGUE) Stellar Param eter Pipeline ()Lee et al.l 
l2007HAllende Prieto et~ani2006[ ). Future versions of the 
pipeline will even estimate [a/Fe]. The unique feature 
of our work is the use of only red and far-red spectral 
regions (6300 A < A < 9100 A). We have chosen this 



spectral region to make use of existing spectra of stars 
in the dwarf satellites of the MW and in the stellar halo 
and dwarf satellites of M31. 

As a general cautionary note on spectroscopic abun- 
dance estimates, studies near the tip of the RGB are 
sensitive to modeling assumptions, including local ther- 
modynamic equilibrium, atmospheric geometry, and the 
mode of energy transport. However, these assumptions 
are significantly less severe than those involved in pho- 
tometric metallicity measurements, and they become 
rapidly less significant at lower luminosities. 

To assess the accuracy and precision of the MRS tech- 
nique, we compare our results to measurements of the 
same stars in the GCs observed with HRS. We discuss 
our observations in §|2|and the preparation of the spectra 
for abundance analysis in §[3l The abundance measure- 
ments are described in §21 We compare the MRS and 
HRS results in §|5l and we quantify errors in [Fe/H]j^j^g 
and [a/Fe]j^^j^g in §[6l We discuss the range of applica- 
tions for the MRS technique, the ways it can improve, 
and its potential advantages over other medium resolu- 
tion techniques in §|7| 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We present observations of seven GCs with the 
Deep Imaging M ulti-Object Spectrometer (deimos, 
iFaber et alll2003[ ) on the Keck ii telescope. Table [T] lists 
the clusters observed, the dates of observations, the ex- 
posure times, and the number of stars targeted. Because 
the targets are very bright, we observed near 12° dusk 
twilight without focusing the primary mirror. Subse- 
quent focusing resulted in insignificant mirror segment 
alignment. The seeing was very good, near 0.5" for 
most observations. We used the OG550 order-blocking 
filter with the 1200 lines mm~^ grating and 0.7" slit 
widths. This configuratio n mimics the DEIMOS confi gura- 
tions of M31 RGB stars (jGuhathakurta et al.ll2006l ) and 
RGB stars in Leo I, a remote dSph satellite of the MW 
(Soh n et aDl2007| ). The spectral resolution is -- 1.3 A 
FWHM (resolving power R ~ 6000). The spectral range 
is about 6300-9100 A with variation depending on the 
slit's location along the dispersion axis. Exposures of 
Kr, Ne, Ar, and Xe arc lamps provided wavelength cal- 
ibration, and exposures of a quartz lamp provided flat 
fielding. 

One 16' X 4' slitmask was devoted to each cluster. Each 
slit included one star. We attempted to maximize the 
number of target stars previously observed with HRS. 
(The slitmask for M79 is an exception, and it contains 
no targets previously observed with HRS.) We selected 
the remaining targets based on CMDs. In order of prior- 
ity, we filled each slitmask with stars from the (1) upper 
RGB, (2) lower RGB, (3) red clump, and (4) blue HB. 
We filled the slits at the edges of the slitmask far from 
the center of the GCs with objects with similar colors 
and magnitudes to stars on the RGB. 

2.2. Data reduction 

We reduce the raw frames using version 1.1.4 of 
the DEIMOS data reduction pipeli ne developed by the 
deep2 Galaxy Redshift Survey.^ IGuhathakurta et "aTl 

^ |http : //astro ■ berkeley . edu/$" {\slm} $cooper/deep/spec2d/ 1 
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TABLE 1 








Observations 






Object 


Date 


Airmass 


Exposures 


# targets 


NGC 1904 (M79)^ 


2006 Feb 2 


1.42 


2 X 300 s 


22 


NGC 2419'' 


2006 Feb 2 


1.21 


4 X 300 s 


70 


NGC 6205 {M13) 


2007 Oct 12 


1.35 


3 X 300 s 


93 


NGC 6838 {M71) 


2007 Nov 13 


1.09 


3 X 300 s 


104 


NGC 7006 


2007 Nov 15 


1.01 


2 X 300 s 


105 


NGC 7078 (MIS) 


2007 Nov 14 


1.01 


2 X 300 s 


63 


NGC 7492 


2007 Nov 15 


1.30 


2 X 210 s 


38 



' ISimon &: Gehal 1)20071 ) have generously provided these observations. 



(|2006f ) give the details of the data reduction. We also 
make use of th e optim izations to the code described in 
ISimon fc Gehal ()2007l . § 2.2 of their article). These mod- 
ifications provide better extraction of unresolved stellar 
sources. 

In summary, the pipeline traces the edges of slits in the 
flat field to determine the CCD location of each slit. The 
wavelength solution is given by a polynomial fit to the 
CCD pixel locations of arc lamp lines. Each exposure of 
stellar targets is rectified and then sky-subtracted based 
on a B-spline model of the night sky emission lines. Next, 
the exposures are combined with cosmic ray rejection 
into one two-dimensional spectrum. Finally, the one- 
dimensional stellar spectrum is extracted from a small 
spatial window in the two-dimensional spectrum encom- 
passing the light of the star. The product of the pipeline 
is a wavelength-calibrated, sky-subtracted, cosmic ray- 
cleaned, one-dimensional spectrum for each target. 

3. SPECTRAL ANALYSIS 

3.1. Determination of spectral resolution 

When we compare model spectra to observed spec- 
tra, we must match the synthetic and observed resolving 
power R. For our configuration of deimos, R is a slight 
function of wavelength, typically varying from 5,500 to 
7,200. We determine the observed spectrum's resolution 
by fitting Gaussians to hundreds of sky lines in the same 
slit as the object spectrum. Then, we fit a parabola to 
the Gaussian widths as a function of wavelength. For 
some short slits, this procedure fails. In those cases, we 
simply fit a parabola to the measured Gaussian widths 
of sky lines from all slits on the same DEIMOS slitmask 
as a function of observed wavelength. 

3.2. Telluric absorption correction 

The last step in obtaining a stellar-only spectrum is 
removal of terrestrial atmospheric absorption. We build 
a telluric absorption template from a continuum-divided 
spectrum of a hot star free of metal absorption lines. On 
2007 November 14, we observed the white dwarf spec- 
trophotometric standard BD-|-28° 4211 with a longslit in 
the same spectrometric configuration as the slitmasks. 
The airmass was 1.018. We assume that all detectable 
absorption lines in the spectrum except Ha are telluric. 

The spectral regions most susceptible to telluric ab- 
sorption are 6864 A < A < 7020 A (B band), 7162 A < 
A < 7350 A, 7591 A < A < 7703 A (A band), 8128 A < 
A < 8352 A, and A > 8938 A. In order to normalize the 
continuum of the telluric absorption template, we simply 
fit a line to the 100 A bands on either side of each region. 



Then, we divide each region by its best-fit line. Because 
the hot star shows no detectable telluric absorption out- 
side of these regions, we set all remaining pixels to 1 to 
prevent introducing noise during telluric absorption re- 
moval. 

We interpolate the telluric absorption template onto 
the wavelength array for each star observed. Then, we 
divide the observed spectrum by the template adjusted 
by the ratio of the airmasses, following the Beer-Lambert 
law: 



(1) 



^Jf(obs)/X(tell) 



where X is the airmass, d is the telluric-corrected spec- 
trum, s is the raw spectrum, and t is the telluric absorp- 
tion template. We carefully adjust the observed spec- 
trum's variance array, treating the noise in the telluric 
spectrum as uncorrelated with the noise in the observed 
spectrum. 

The amount of attenuation varies as a function of 
the number density of absorbers — particularly water 
vapor — along the line of sight. The associated frequency- 
dependent optical depth varies on timescales possibly less 
than an hour. Regardless, the telluric absorption removal 
procedure described here works very well even for ob- 
servations taken 21 months before the telluric standard. 
The procedure does a poor job in three spectral regions: 
the saturated A band, the saturated portion of the B 
band, and a small 40 A region around 8245 A. We ignore 
these regions in the abundance analysis. 

Figure [1] demonstrates the efficacy of the telluric ab- 
sorption removal. The upper panel shows part of an ex- 
ample spectrum before the telluric correction, along with 
the telluric absorption template. The middle panel shows 
the example spectrum after the correction. The bottom 
panel includes a synthetic stellar spectrum (see details 
in §[3]), which shows that every absorption line left after 
the telluric correction is intrinsic to the star. 

3.3. Radial velocities 

We measure radial velocities of each target star to 
check its cluster membership and to shift its spectrum 
to the rest frame for abundance measurement. We cross- 
correlate each telluric-divided observed spectrum with a 
synthetic spectrum (see details in §|4l) with Toff — 4500 K, 
\ogg = 1.5, [Fe/H] = -1.5, and [a/Fe] = -hO.2 in the 
spectral region 8450 A < A < 8700 A. This region in- 
cludes the Ca II triplet, which is strong even in hot, ex- 
tremely metal-poor stars, making it ideal for radial ve- 
locity determination in a wide range of stars. We shift 
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Fig. 1. — A section of the spectrum for one star in M13. Top: 
The observed spectrum (black) with the telluric absorption tem- 
plate (green). Every feature in the telluric absorption template has 
a counterpart in the observed spectrum, but the reverse is not true. 
Middle: Applying the telluric absorption correction completely re- 
moves some of the absorption features seen in the raw spectrum, 
makes a few of the absorption features weaker, and leaves many 
of the features unaffected. Bottom: The telluric-divided spectrum 
(black, same as in middle panel) with the best-fit synthetic spec- 
trum (red). The close agreement between the two spectra demon- 
strates that the features remaining in the observed spectrum are 
intrinsic to the star. Both the green and red lines have been scaled 
by the measured continuum. 

the observed spectra to the rest frame to complete the 
remainder of the analysis. 

3.4. Continuum determination 

The abundance measurements are particularly sensi- 
tive to an accurate determination of the continuum. Un- 
derestimating the continuum will make absorption lines 
appear too shallow, and the derived abundances will be 
too low, the temperature too high, or both. The abun- 
dance analysis described in §|4]is insensitive to the global 
continuum shape and instead relies on high-frequency 
line-to-line variations in flux as a function of wavelength. 
Therefore, we have decided to focus on local continuum 
determination. 

First, we determine the spectral regions free of stellar 
absorption. Following the procedures in §|4l we generate 
a synthetic spectrum between 6300 A and 9100 A of a 
star with T^s = 4300 K, logg = 1.5, and [Fe/H] = -0.5. 
The spectrum is smoothed through a moving Gaussian 
kernel of tr = 0.6 A to simulate the approximate spec- 
tral resolution of deimos. The synthetic spectrum has 
a perfectly flat continuum, and the units are such that 
the continuum is 1. We call spectral regions with syn- 
thetic flux greater than 0.96 and a minimum width of 
0.5 A "continuum regions." Pixels at observed wave- 
lengths outside of these regions will not contribute to 
the continuum determination of observed spectra. 

Next, we compute the continuum. Each pixel in the 
continuum (cj ) is the weighted average of its neighboring 
pixels in the observed spectrum (si). The weight is a 
combination of the inverse variance {af) and proximity 
in wavelength. 
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Fig. 2. — A different section of the same star as in Fig. [T] 
Top: The observed spectrum (black) and the measured contin- 
uum (green). Middle: The observed spectrum divided by the con- 
tinuum. Bottom: The continuum-normalized observed spectrum 
(black) with the best-fit synthetic spectrum (red) , wh ich h as a lso 
been subjected to continuum normalization (see ^ 13.41 and S I4.4|I . 



1 if Xi G continuum region 
otherwise 



(4) 



This process simultaneously accomplishes smoothing and 
interpolation across non-continuum regions. We normal- 
ize each observed spectrum by dividing by its measured 
continuum. Figure [2] shows the measured continuum in 
one example spectrum. 

The continuum regions include some weak absorption 
lines, which will drag the continuum determination down 
by up to 4% in the coolest, most metal rich stars in 
our sample. In S I4.4( we determine the synthetic spec- 
trum that most closely matches the observed spectrum. 
Before comparing a synthetic spectrum to the observed 
one, we apply the above continuum determination tech- 
nique to both. Even though the synthetic spectrum has 
a perfectly flat continuum, application of this continuum 
determination technique causes the same degree of weak 
feature continuum suppression in both the synthetic and 
observed spectra, thereby correcting this effect. 

3.5. Pixel mask 

Each spectrum has its own pixel mask. Table [2] lists 
spectral regions in observed and rest wavelength that 
may confuse the abundance analysis. The first part of 
the table lists regions of saturated or very strong telluric 
absorption. We mask the pixels that fall in these regions 
before we shift the spectrum to zero velocity. The second 
part of the table lists stellar absorption features that are 
difficult to model. We mask the pixels that fall in these 
regions after shifting the spectrum to zero velocity. These 
features include the strongest a element lines in the spec- 
tral range of deimoS: the Ca II triplet and Mg I A8807. 
The cores of the Ca II lines are formed very high in the 
stellar atmosphere out of local thermodynamic equilib- 
rium (LTE). Because we could not synthetically repro- 
duce the width and shape of the Mg I line in Arcturus 
(see § l4.2p . we excluded it from the abundance analysis. 

In addition to the standard mask in Table [21 we in- 
spected each spectrum to mask improperly subtracted 
sky lines, cosmic rays, and other obvious instrumental 
artifacts. 

4. CHEMICAL ABUNDANCE ANALYSIS 
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TABLE 2 
Spectral Masks 



Feature 


Lower Wavelength 


(A) 


Upper Wavelength (A) 




Telluric Absorption 


B band 
A band 
strong telluric absorption 


6864 
7591 
8225 




7020 
7703 
8265 


Stellar Absorption 



Ca I A6343 


6341 


6346 


Ca I A6362 


6356 


6365 


Ha 


6559.797 


6565.797 


K I A7665 


7662 


7668 


V I A8116,A8119 hyperfine structure 


8113 


8123 


poorly modeled absorption in Arcturus 


8317 


8330 


Ca II A8498 


8488.023 


8508.023 


Ca II A8542 


8525.091 


8561.091 


Ca II A8662 


8645.141 


8679.141 


Mg I A8807 


8804.756 


8809.756 



TABLE 3 
Grid of Synthetic Spectra 



Parameter Lower Limit 



Upper Limit 



Step 



Number 



ToH (K) 



4000 





( 0.0 


4.5 




0.0 


5.0 




0.5 


5.0 




0.5 


5.0 


log 9 S 


1.0 


5.0 




0.5 


5.0 




1.0 


5.0 




1.0 


5.0 




I 1.5 


5.0 


Fe/H] 


-4.0 




a/Fc] 


-0.6 





8000 

f Toff < 4500 and [Fe/H] < -2.5 

f (Tcft < 4500 and [Fe/H] > -2.5) or 4500 < T^s < 6000 
f 6000 < Tcff < 7000 
f 7000 < Tcff < 7200 and [Fe/H] > -3.5 
f 7000 < T^ff < 7200 and [Fe/H] < -3.5 
f 7200 < Tcft < 7500 and [Fe/H] > -2.5 
f 7200 < Tcff < 7500 and [Fe/H] < -2.5 
f Teff > 7500 and [Fe/H] > -2.5 
f Teff > 7500 and [Fe/H] < -2.5 
0.0 
+1.0 



100 if Tcff < 5500 
200 if %ft > 5600 



29 





r 10 




11 




10 




10 


0.5 I 


9 




10 




9 




9 




I 8 


0.1 


41 


0.1 


17 



Total 



210 902 



The heart of this analysis is a very large grid of syn- 
thetic stellar spectra. The grid contains four dimen- 
sions: effective temperature (Tofi), surface gravity (logg), 
metallicity ([Fe/H]), and alpha enhancement ([a/Fe]). 
Table [3] lists the details of the grid. 

In this article, we use the standard spectroscopic abun- 
dance notation. The ratio of any two elements A and B 
in a star relative to their ratio in the Sun is 

[A/B] = log[n(A)/n(B)] - log[n(A)/n(B)]0 (5) 

where n is number density.^ [Fe/H] represents metal- 
licity, the overall bulk content of the heavy elements. 
Because neutral Fe and neutral a element lines domi- 
nate the red spectra of stars in our reff-log(7-[Fe/H] do- 
main, we assume [X/H] ~ [Fe/H] for most elements. We 
stress, however, that lines of all elements except the a 
elements contribute to the MRS measurement of [Fe/H] . 
The abundances of the a elements Mg, Si, S, Ar, Ca, and 
Ti are modified by the additional parameter [a/Fe] : 

^ In this article 12 + log[n(Fe)/n{H)]0 = 7.52 (as adopted by 
ISneden et al.in"992l). The abundan ces of all other elements are the 
solar values of lAnders &: Grevess3 1)19891 ) except Li, Be, and B, for 
which we use the meteoritic values. 



[a/H] = [Fe/H] + [a/Fe] 



(6) 



The abundance enhancements of all six elements vary 
together. 

We have found the best agreement with HRS by fix- 
ing the microturbulent velocity (vt) to logg with an em- 
pirical relation. The best linea r fit to the RGB sample 
(log 5 < 3.3) of lFulbrightJ HQOS) is 



Vt = 2.700 - 0.509 logg km/s. 
4.1. Model atmospheres 



(7) 



The starting point for generating a synthetic spec- 
trum is a model stellar atmosphere, which is a tabu- 
lation of temperature, pressure, electron fraction, and 
opacity as a function of optical depth. We choose to use 
the Castelli fc Kurucz grid of mo dels with no convective 
overshooting (jCastelli et aLlflOOTI ) . For atmospheres with 
[Fe/H] > -2.5, we use the "ODFNE W" models, with up- 
dated opacity distribution functions (|Castelli fc Kuruca 
[200l . 

The published model atmosphere grid points of T^s, 
[Fe/H], and [a/Fe] are coarser than the step sizes listed 
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in Table [3l Therefore, for a fixed array of optical depths, 
we linearly interpolate temperatures, pressures, electron 
fractions, and opacities between model atmosphere grid 
points to generate a single model atmosphere for every 
grid point in Table [3l 

The alpha enhancement changes the ionization balance 
and free election fraction in the atmosphere. ODFNEW 
([Fe/H] > -2.5) models are available for [a/Fe] = 0.0 
and +0.4. For synthesizing spectra with [a/Fe] < 0.0, we 
choose the former and for spectra with [a/Fe] > +0.4, 
we choose the latter. For intermediate values of [a/Fe], 
we linearly interpolate between the two regimes. At- 
mospheres with [Fe/H] < —2.5 are available only with 
[a/Fe] — 0.0. For the model atmospheres only, the a 
elements are O, Ne, Mg, Si, S, Ar, Ca, and Ti. 



4.2. Line lists 

We assembled a line list of wavelengths, excitation po- 
tentials (EP), and oscillator strengths {gf) for atomic 
and molecular transitions that occur in the red spectral 
regions of stars in our stellar parameter range. In stan- 
dard HRS analyses, it is often possible to use only transi- 
tions with accurately measured oscillator strengths. For 
this MRS analysis that covers a broad spectral range with 
many blended lines, we must use a database of lines with 
oscillator strengths of varying accuracy. 

To beg in, we queried the V ienna Atomic Line Database 
(VALD, iKupka et al.l Il999f ) for all transitions of neu- 
tral or singly ionized atoms with EP < 10 eV and 
\oggf > —5. We supplemen ted the l ist with CN, C2, and 
MgH molecular transitions (jKurucdll992) and Li, Sc, V , 
Mn, Co, Cu, and Eu hyperfine transitions ()Kurucdll993l ). 
By far, CN is the most important molecule for our red 
spectra. TiO is also strongly present in cool, metal-rich 
stars. However, the TiO system is extremely complex 
and difhcult to model accurately. The large number of 
TiO electronic transitions in red spectra makes TiO spec- 
tral synthesis computationally daunting. Fortunately, 
TiO becomes a significant absorber only in metal-rich 
stars with Tcff < 4000 K. Therefore, we did not include 
TiO in our line list, and we chose Tcff = 4000 K as the 
grid's lower limit. 

Next, we generated synthetic spectra of the Sun 
and Arcturus using our line list, model atmospheres 
as described in ij |4.1i and the current ver sion of 
the LTE spectral synthesis software MOOG IjSnedenI 
|1973[) . The spectral range is 6300 A to 9100 A, 
and the resolution is 0.02 A. The line broadening ac- 
counts for collisions with neutral h ydrog en for for lines 
with tabulated damping constaiits (IBarklem et al.ll2000l : 
iBarklem fc AsDelund-Johansson|[2005D . For other lines, 
the Unsold approximation multiplied by 6.3 gives the van 
der Waals line damping parameter, and MOOG calcu- 
lates additional radiative and Stark broadening. How- 
ever, the choice of damping parameter does not notice- 
ably affect the spectra of low-pressure giant stars such 
as Arcturus, nor any star in the domain of Tcs and \ogg 
that we consider here. 

The atmospheric parameters for the Sun are Tcff = 
5798 K, log 5 = 4.44, and [Fe/H] = 0.0. For Arcturus, we 
adopt the atmospheric pa rameters and non-sola r abun- 
dance ratios determined bv lPeterson et al.1 (|1993D : Tcff = 



TABLE 4 
Line List 



Wavelength (A) 


Species 


EP (eV) 


logs/ 


8647.703 


22.0 


4.654 


-3.561 


8647.792 


25.1 


6.834 


-2.489 


8647.799 


21.0 


4.558 


-1.759 


8647.807 


26.0 


5.720 


-1.519 


8647.940 


26.0 


6.5 


-0.5 


8648.024 


21.0 


4.126 


-2.710 


8648.040 


607.0 


1.093 


-3.058 


8648.366 


607.0 


1.093 


-1.096 


8648.400 


607.0 


0.912 


-2.763 


8648.455 


14.0 


6.21 


0.0 


8648.556 


16.0 


8.408 


-1.720 


8648.630 


27.059 


2.280 


-3.852 


8648.688 


607.0 


0.912 


-1.466 


8648.699 


27.059 


2.280 


-4.068 


8648.716 


20.0 


4.554 


-2.024 



Note. — Table |4] is published in its entirety in 
the electronic edition of the Astrophysical Jour- 
nal. A portion is shown here for guidance regard- 
ing form and content. The first column is the line 
wavelength. The second is a code representing the 
atomic or molecular species. The integer part of 
the code is the element's atomic number. Codes 
greater than 100 represent molecules. For exam- 
ple, 607 represents CN. The first decimal place of 
the code is the ionization state, where is neutral. 
Remaining decimal places are mass numbers for 
isotopic hyperfine transitions. The third column 
is EP. The invented Fe I transitions have EP with 
fewer than three decimal places. The fourth col- 
umn is oscillator strength. The transitions that 
we modified have oscillator strengths with fewer 
than three decimal places. 

4300 K, logg = 1.50, and [Fe/H] = -0.5 0.^ We com- 
ared these two synthetic spectra to the iHinkle et al.l 
2000( ) atlases of the Sun and Arcturus. After smoothing 
the synthetic spectra through a Gaussian kernel to match 
the line profiles of the atlas, we inspected the spectra in 
detail. First, we adjusted oscillator strengths of aber- 
rant atomic lines in the solar synthetic spectrum until 
the strengths of the synthesized lines matched that of 
the observed lines. Then, we repeated the process for 
Arcturus, making adjustments that did not cause dis- 
agreement in the solar spectrum. 

Occasionally, we encountered observed lines not 
present in the line list. Most of these cases could be 
resolved by relaxing the log gf > —5 restriction. How- 
ever, in cases where a single line was not represented in 
the VALD or Kurucz line lists, we invented a transition 
of Fe I of the EP and gf required to reproduce the line's 
observed strength in both the Sun and Arcturus. The 
final line list contains 30873 atomic and 17345 molecular 
transitions, representing 71 elements. Tabic 0] presents 
the entire line list. 

With the final line list, the mean absolute deviation 
between the pixels of the solar spectrum and the pixels 
of its synthesis is 5.6 x 10~^, and the standard deviation 
is 1.3 X 10^^. For Arcturus, the mean absolute deviation 



^ I Peterson et al.l II1993I ) use a value of the solar iron abundance 
that is 0.15 dex higher than the one adopted in this article. In 
fact, they note that decreasing the iron abundance to the value 
used here reproduces the profiles of weak Fe I lines without adjust- 
ing laboratory-measured oscillator strengths. We find that their 
published value of [Fe/H] with our value of the solar iron abun- 
dance reproduces the spectrum of Arcturus very well. 



Abundances from Medium Resolution Spectroscopy 



is 9.4 X 10 ^, and the standard deviation is 2.1 x 10 ^. 
The units are such that the continuum is 1. 

4.3. Synthetic spectra generation 

We generate the hbrary of synthetic spectra from the 
final line list. Each spectrum ranges from 6300 A to 
9100 A with a resolution of 0.02 A. To save computa- 
tion time in the abundance determination, we bin each 
synthetic spectrum by a factor of 7 (0.14 A resolution). 
When the binned and unbinned spectra are smoothed 
to the best resolution provided by DEIMOS (~ 1.1 A 
FWHM), individual pixels in the binned spectra differ 
from the unbinned spectra by less than one part in 10'^. 

4.4. [Fe/H] and [a/FeJ determination 

Our continuum normalization for observed spectra (see 
i) l3.4|) may be depressed very slightly by weak absorption 
lines. To counteract this effect, we allow weak absorp- 
tion in synthetic spectra to affect their normalizations in 
the same way that they affect observed spectra. First, 
we interpolate a synthetic spectrum onto the same wave- 
length array as the observed spectrum. Then, we smooth 
the synthetic spectrum through a Gaussian filter whose 
width is the observed spectrum's measured spectral res- 
olution as a function of wavelength (see § 12. 2p . To com- 
plete the renormalization, we divide the synthetic spec- 
trum by the synthetic "continuum" determined exactly 
the same way as in § 13.41 including the continuum regions 
and weighting by the inverse variance of the observed 
spectrum. 

We discard all stars with photometric Tcs < 4000 K 
because such stars are susceptible to TiO absorption, 
which we do not attempt to model. The coolest tem- 
perature on the synthetic grid is 4000 K for this reason. 

Next, we compute x^ for an initial guess at the four at- 
mospheric parameters T^s, logg, [Fe/H], and [a/Fe]. A 
Levenberg-Marquardt optimization algorithm finds the 
values within the bounds of the grid that minimize x^ 
computed from the difference between the observed spec- 
trum and trial synthetic spectrum. We sample the pa- 
rameter space between grid points by linearly interpolat- 
ing the synthetic spectra at the neighboring grid points. 
The parameters for the best-fit spectrum are the final 
atmospheric parameters for the star. The fitting error 
on each parameter is \/ o'^y^ where a^ is the parameter's 
diagonal element of the covariance matrix. The fitting 
errors are usually much smaller than the total error, in- 
cluding systematic error, which will be discussed in §[S1 
Fitting errors are particularly small for high S/N spectra. 

Finding the global minimum x^ in four dimensions is 
difficult. Surface gravity is particularly difficult to mea- 
sure. First, the red spectral region contains very few 
gravity sensitive lines that are easy to model. Second, 
few transitions from ionized species are visible in our 
spectra, making surface gravity poorly constrained. To 
reduce the dimensionality of the parameter space, we fix 
Teff and logg based on photometry (P. B. Stetson, pri- 
vate communication) and theoretical isochrones shifted 
to the distance modulus of the target. MRS abundances 
almost always agree with HRS abundance s more c losely 
when we fix Teff and log g photometrically. [Harris! ()1996L 
revised 2003) provides the distance modulus and redden- 
ing of each cluster. 



For the GCs in this article, we use 14.0 Gyr, 
[g/Fe ] = -1-0.3 isochrones f rom the Yonsei-Yale group 
(YY, iDemargue et al.l I2004D . Teg and log g are mostly 
insensitive to age and ot enhancement. Given a star's 
/ magnitude and V — 1 color, we linearly interpo- 
late between tracks of constant age to recover Toff and 
log 5. (We do not have / measurements for a few 
stars in M13 nor any star in NGC 7492. Instead, we 
use V and B — ]/.) We have also exper imented with 
Victoria-Regina ([VandenBerg et all [2006) and Padova 
(jGirardi et al.l 12002) isochrones. Although the Padova 
photometric [Fe/H] is quite different from the other two 
sets of isochrones near the tip of the RGB, the pho- 
tometric Teff and logg do not change enough to affect 
abundance analysis. We have also compared theoreti- 
cal YY temperatures to th e empirical temperatures of 
iRamirez fc MelendezI ()20Q5l . RM05). For Teff < 4800 K, 
the average discrepancy AToff < 50 K. For T^s > 5000 K, 
the RM05 Tes tends to be up to 250 K lower than the YY 
Toff. Given their success in reproducing GC abundances 
(§ 15. 2p . we have chosen to use YY temperatures. 

Future applications of this method to inhomogeneous 
stellar populations may require an assumption of the 
age or iteration of the age until the spectroscopic [Fe/H] 
agrees with the photometric [Fe/H]. In theory, this iter- 
ation will give the age of the star. However, photometric 
errors and the systematic errors in both abundance mea- 
surement methods will undoubtedly make such an age 
very uncertain. In reality, Teff and logg change little 
with age or metallicity for a given V — I color, so the as- 
sumed isochrone parameters hardly affect the abundance 
results. 

5. RESULTS AND COMPARISON TO HRS 
ABUNDANCES 

To demonstrate the effectiveness of the MRS abun- 
dance technique, we have observed GCs specifically to 
compare MRS to HRS metallicities and a enhancements. 
In this section, we present those comparisons and exam- 
ine MRS results for systematic trends. 

5.1. Toff 

The MRS technique can recover T^g even without pho- 
tometry. Figure [3] shows the Toff of the best-fit atmo- 
spheric parameters when we allow Toff, log g, [Fe/H], and 
[a/Fe] to vary compared to the published Toff used in 
HRS studies of the same stars. Not all authors of HRS 
studies determine Toff spectroscopically. The figure leg- 
end indicates which HRS studies choose T^s such that 
lines of different excitation potential yield the same abun- 
dance ("spec") and which HRS studies rely exclusively 
on broadband photometry to determine Toff ( "phot" ) . 

(Tcff)MRS reproduces (Toff)HRS very well in most cases. 
The absolute deviation for 77% of stars falls within 
150 K. At Toff > 5000 K, MRS analysis underpredicts 
Toff because at high temperatures, lines with low EPs be- 
come immeasurably weak. [An alter native explanation is 
that th e photometric temperatures of lCohen fc Melended 
()2005a| ) could be inaccurate.] Spectroscopic temperature 
is measured by comparing the strengths of lines with a 
range of EP. Without the low-EP lines, the temperature 
is difficult to measure. Fixing Toff photometrically alle- 
viates this problem and also eliminates the large random 
error in a spectroscopic Toff of lower S/N stars. For the 
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Fig. 3. — T^g determined from MRS without photometric infor- 
mation for individual stars (this work) vs. T^ff used in HRS studies 
of the same stars. (T'j.ff )hrs is not necessarily determined spectro- 
scopically. The figure legend indicates whether (rcff)HRS is de- 
termined spectroscopically ("spec") or photometrically ("phot"). 
At T(,ff > 5000 K, MRS analysis may underpredict Tog be- 
cause at high temperatures, lines with low EPs become immea- 
surably weak. Alternatively, the photometric temperatures for 
M13 based on theoretical atmospheres (Cohen & Melcndez 2005a) 
may be overestimates. The dotted line is one-to-one. The follow- 
ing refer ences provide (Tpff)Hps- SC SOl: Shetronc et al. (2001), 
CM05a: ICohen & Mel^nded lJ2005a|), S04: Sncdcn et al. (2004), 
RC02: Ramirez & Cohen (2002), MPS03: Mishenina et al. (2003), 
K98: Kraft et al. ( 1998) , S97: Sneden et al. (1997). SPK OO: 
ISneden et al.l 1 12000) . and CM05b: ICohen fc MelendezI (I2005bl) . 
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Fig. 4. 



[Fe/H], 



Teff determined from V — I color of 



individual stars in Ml3 and M15. (A few stars in M13 rely on 
B — V rather than V — I.) The lack of a trend suggests little 
systematic covariance between measur ements of [F e/Hl^^rpg and 
Tcft. The dotted lines show [Fe/H]jjj^g UPritzl et al.ll2005D . 



remainder of this article, both Teg and log 5 are set by 
photometry. 

5.2. [Fe/H] 

We expect any spread in [Fe/H] within a single GC to 
be random error because these GCs are monometallic. 
Correlation of [Fe/H] with other parameters can indicate 
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Fig. 5. — [F^/H]iy[jj3 vs. [Q!/Fe]]y[jj^g of individual stars in M13 
and M15. The lack of a trend suggests little systematic covariance 
between measurements of [g/Fej^ ^nc; vs . [Fe/H] j^j j^^g . The dotted 



between measurements of fo/Fe 
lines show [Fe/H]jjp^g (.Pritzl et 



alj |2005l) . 



systematic errors. First, we show in Fig. [¥]the spectro- 
scopic [Fe/H] vs. photometric Toff in M13 and M15. In 
all of M15 and above 5000 K in M13, we see only random 
scatter. The bright, cool stars near the tip of the RGB in 
M13 show a slight positive slope of ~ 0.1 dex from 4000 K 
to 5000 K. The covariance is undesirable but not unex- 
pected. Both higher temperatures and lower metallicities 
weaken absorption lines. The magnitude and significance 
of the trend are small and restricted to the hottest stars. 
Furthermore, we emphasize that M13 is the worst case. 
No other GC shows this correlation. 

We also expect no correlation between [a/Fe] and 
[Fe/H] within a single GC. Both [a/Fe] and [Fe/H] affect 
the strength of a element absorption lines. If there were 
a systematic trend, we would expect the two parameters 
to be anti-correlated. No cluster shows a convincing anti- 
correlation. Figure [5] shows both parameters for M13 and 
M15. M13 is again the worst case, and the significance of 
any trend is destroyed by removing the one or two points 
with the most discrepant [Fe/H]. 

The MRS technique matches the metallicities deter- 
mined from HRS to within the HRS metallicity scat- 
ter for a given GC. Figure [6] shows [Fe/H] of the best- 
fit atmospheric parameters (where photometry deter- 
mines Teff and \ogg) compared to [Fe/H]jjj^g determined 
from Fe I lines in the range -2.38 < [Fe/H] < -0.76. 
We choose Fe I because Fe I lines greatly outnumber 
Fe II lines in the deimos spectra with our range of 
Tcff and logg. Although these GCs are almost certainly 
monometallic, we have plotted the measurements for in- 
dividual stars on both axes. The range of [Fe/H]jjj^g 
within a single GC, particularly depending on the au- 
thors of the measurement, demonstrates the uncertainty 
to which the mean metallicity is known. The MRS mea- 
surements do not fall outside of this scatter. 

We measure individual metallicities for all stars that 
we observe. We eliminate non-members by radial veloc- 
ity, and we discard stars with photometric Tcff outside 
the range of the spectral grid. We also discard HB stars 
and stars with noticeable TiO absorption. Figure[7]shows 
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TABLE 5 
Average Globular Cluster Abundances Comparison 



Cluster 


A^HRS^' 


[Fe/Hlpvios'^ 


(["/Fc]hrs>" 


A^DEIMOS'* 


([Fe/H]MRs> 


rms([Fc/H]MRs) 


(["/FcImrs) 


rms([a/Fe]njR3) 


M79 


2 


-1.42 


+0.26 ± 0.08 


33 


-1.69 + 0.02 


0.15 


+0.27 + 0.01 


0.11 


NGC 2419 


1 


-2.32 


+0.20 ±0.08 


30 


-2.28 + 0.02 


0.08 


+0.28 + 0.01 


0.09 


M13 


60 


-1.57 


+0.20 + 0.01 


69 


-1.66 + 0.02 


0.11 


+0.25 + 0.01 


0.09 


M71 


40 


-0.76 


+0.27 + 0.01 


47 


-0.92 + 0.04 


0.17 


+0.27 + 0.01 


0.17 


NGC 7006 


6 


-1.55 


+0.24 + 0.01 


20 


-1.59 + 0.03 


0.13 


+0.35 + 0.01 


0.08 


M15 


49 


-2.38 


+0.42 + 0.03 <= 


44 


-2.42 + 0.01 


0.12 


+0.33 + 0.01 


0.07 


NGC 7492 


4 


-1.85 


+0.23 + 0.02 


21 


-2.04 + 0.02 


0.08 


+0.32 + 0.02 


0.10 



Note. — The rms values and errors on the mean are weighted by the inverse square of individual measurement errors. 
^ Number of stars observed with high resolution spectroscopy.'' Pritzl ct al. (2003), high resolution spectroscopy.*^ See Fig. [7]for references 
from which these averages were calculated. "* Number of DEIMOS spectra analyzed in this article. ° See i) l5.3l for a discussion of the anomalously 
large value of [o/Fe]jjp^g for M15. 
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Fig. 7. — Mean [Fe/H] values in GCs observed with deim os vs. 
values of [Fe/H] determined from HRS (|Pritzlet al.l[2005t ). The 
dotted line is one-to-one. The full error bars are weighted sam- 
ple standard deviations, and the hash marks along the error bars 
are weighted errors on the mean. The HRS errors are weighted 
by the inverse square of the measurement err ors of all of the indi- 
vidual stars in the following references. M79: IGratton &: Ortolanil 
(1989), NGC 2419 (the HRS error bar is the individual measure - 
ment error for the only star available): [S hetrone et al. (2001|), 
M13: ISneden et al.1 lf20M') andlCohen fc Mel cndcz (2005a), M71: 
IGratton et al .r (19861. ISnedenet al.l l|199l) . IRamire z & Cohe3 
(|2002l) . and" Mishenin a et al.l 11200 J) . NG C 7006: [K raft et ^ 
jTsm ). M15: jSneden et al. fl997) and ISneden et al.l (l2000ir 



Fig. 6. — This work's [Fe/H] vs. [Fe I/H] determined from neutral 
iron lines in HRS. The dotted line is one-to-one. The symbols, 
colors, and references are the same as in Fig. [3] 



the weighted mean of [Fe/H]j^j^g of all the stars in each 
cluster . The abscissa is [Fe/H]jjj^g from the compila- 
tion of iPritzl et al.l ()2005l ). Two types of errors bars are 
shown on both axes: the standard error on the mean and 
the sample standard deviation. The latter represents the 
typical measurement error on one star. The HRS error 

bars are determined from [Fc/H]jjj^g of individual stars ^GC 7492: [Cohen fc Melcndcz (2005bJ 

in the references given in the figure caption. The error 
bars for both MRS and HRS are weighted by the inverse 
square of the individual measurement errors. Table [5] 
lists the same data along with the number of stars in 
each sample. The HRS measurements for NGC 2419, 
M79, and NGC 7492 are based respectively on only 1, 2, 
and 4 stars. We suggest that, in these cases, a cluster's 
mean [Fe/H] may be more precisely determined with our 
larger MRS samples. 

Table [7] lists MRS results for individual stars. Where 
available, all three photometric magnitudes (i?, V, and 
/) are given, but only two determine T^b and log 5. V 
and / are preferred. Tcs and log 5 in this table are the 
values that we have determined photometrically. The 
last three columns are data from and references to HRS 
studies for the stars in common between the MRS and 
HRS datasets. 

All GCs seem to display internal variations in [C/H] 



and [N/H] (^e.g.. lGrundahlfl999l) . iCohen et all (l2005h 
show ranges of 2-3 dex of [N/H] in five GCs, includ- 
ing M13, M71, and M15. The variations could alter the 
strength of CN absorption. M71 is the only GC in our 
sample metal-rich enough to exhibit strong red CN ab- 
sorption. If star-to-star CN abundance variations affect 
[Fe/H]jy[j^g, then the metallicity determined only from 

the spectral region affected by CN (7850-8400 A) should 
differ from the metallicity determined from the rest of 
the spectrum. We subjected M71 to this test, and [Fe/H] 
measurements from the two cases agree very closely for 
every star. CN abundance variations, if they exist, do 
not appear to contribute to the error on individual stars 
or to the large scatter in [Fe/H]j^j^g for M71. 

5.3. [a/Fe] 
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Fig. 8. — Top: Portions of an observed spectrum of a star in 
M13 (black) and the best fit synthetic spectrum (red), which has 
[o/Fe] = +0.40. Fe fines are labeled in red and a element lines 
are labeled in green. The residual is the difference between the 
observed and synthetic spectra. Bottom: The same observed spec- 
trum (black) and a synthetic spectrum (red) with the same pa- 
rameters as the synthetic spectrum in the top panel except that 
[a/Fe] = 0.0. The observed spectrum is highly inconsistent with 
the solar value of [a/Fe]. 



]\/IRS can give an abundance dimension beyond metal- 
licity. IVEost of the stellar absorption lines visible in 
DEIMOS spectra are from Fe I, but there are a comparable 
number of a element absorption lines. By far, the most 
abundant are Ti, Si, Ca, and JVIg, in order of decreasing 
prevalence. Figure[H]shows two spectral regions with high 
concentrations of lines from all of these elements except 
Ca. The figure also shows two synthetic spectra that are 
identical except for their values of [a/Fe]. The observed 
spectrum is consistent with the synthetic spectrum for 
which [a/Fe] = -1-0.40, but highly inconsistent with the 
synthetic spectrum for which [a/Fe] = 0.0. The dis- 
crepancy demonstrates that MRS can easily distinguish 
between the halo plateau and solar values of [a/Fe]. 

In order to build additional confidence in the abil- 
ity of ]\/IRS to determine [a/Fe], we compare values 
of [a/Fe],^j^g to HRS measurements. Figure [9] shows 
[a/Fe]]y[j^g compared to HRS measurements of [Mg/Fe], 
[Si/Fe], [Ca/Fe], and [Ti/Fe]. AU four plots show weak 
or no correlations. The lackluster agreement is not sur- 
prising because [o^/FeJj^j^g is a weighted combination of 
all four elements. 

To better compare IVIRS and HRS results, we define 
[a/Fe]jjpjg, which is a weighted mean of the four a ele- 
ment ratios. 



[«/Fe]HRs = (7Mg [Mg/Fe] + 473; [Si/Fe] (8) 

+27ca [Ca/Fe] +67Ti [Ti/Fe]) 
X (7Mg + 47si + 27ca + 67Ti)"^ 
S[a/F^Ws = [7Mg (<5[Mg/Fe])' + 731 (4<5[Si/Fe])' (9) 
+7ca(2(5[Ca/Fe])V7Ti(6(5[Ti/Fe])']'/' 
X (7Mg + 47si + 27ca + 67Ti)"^ 

The 7x coefficients are 1 if the element X has a pub- 
lished HRS value and otherwise. J[X/Fe] is the pub- 
lished error. For those studies that quote [X/H] instead 
of [X/Fe], 6[X/Fe] 



weights 1, 4, 2, and 6 mimic the relative presence of 
absorption from the four elements. This combination of 
weights also gives good correlation between [a/Fe]j^j^g 
and [a/Fe]jjf{g, shown in Fig. [TOl These data are also 
listed in Table [71 Although the weights in [a/FeJjjpg 
approximate the weights that the individual elements re- 
ceive in determining [a/Fe]]yjj^g, matching the elemental 
averages exactly is not possible. Therefore, the compar- 
isons between [a/Fe]j^j^g and [a/Fe]jjj^g are approximate 
checks of agreement. 

Three studies deserve partic u lar m ention. The M15 
measurements of ISneden et al.l ()2000( l use Hydra spec- 
tra, which have lower spectral resolution than most HRS 
studies. They do not attempt to measure Mg, and most 
of their stars do not display enough Si or Ti absorption to 
permit accurate measurements of those elemental abun- 
dances. Furthermore, their Ca measureme nts are based 
on onl y one line. The M13 measurements of lSneden et al.l 
(|2004[ ) do not include Si, Ca, or Ti. They do measure 
Mg, but Mg is the least visible a element in the deimos 
spectra. For th ese reasons, we excl ude t he [a/Fe]^p,g 



^((5[X/H])2-h(<5[Fe/H])2. The 



measu rements of lSneden et al.l ()2000( ) and ISneden et al.l 
(|2004f l from Fig.HOjand Tabl e III We also draw at tention 
to the M15 measurements of lSneden et al.l ()1997f ). They 
measure particularly large values for [Si/Fe], and they 
do not recover [Ti/Fe] for any of the stars that are in 
common between their HRS sample and our MRS sam- 
ple. Therefore, [ck/FeJjjpg has a different meaning for the 
M15 data (green points in Fig.[TO| than for the data from 
other clusters. 

Comparing individual stars limits the sample to stars 
observed with both MRS and HRS. In order to draw on 
the full set of stars observed with either MRS or HRS, 
we compute mean values of [a/Fe] for entire clusters. 
([a/Fe]) is weighted by the inverse square of 5 [a/Fe]. 
Figure [TT] shows the comparison. These averages are also 
listed in Tabled The range of ([a/Fe]) for these GCs 
is very smaU, but ([a/Fe])MRS and ([a/Fe])HRS agree to 
< 0.1 dex for all clusters. 

Finally, we show ([a/Fe]) vs. ([Fe/H]) in Fig.[12l Also 
shown is the relationshi p between these tw o quantities for 
MW stars tabulated bv lVenn et all (|2004D . The ([Fe/H]) 
values from both MRS and HRS are consistent with the 
halo plateau value of [a/Fe] ^ -f 0.3. 



6. QUANTIFICATION OF ABUNDANCE ERRORS 
6.1. Total Error on [Fe/H] 

The Levenberg-Marquardt algorithm which determines 
the best-fit synthetic spectrum by minimizing x^ gives 
an estimate of the fitting error based on the depth of 
the x^ minimum in parameter space. However, the fit- 
ting error is usually a small part of the total error. The 
major sources of error in high S/N spectra are errors 
in atmospheric parameters and imperfect spectral mod- 
eling. Two effects make the total error a function of 
[Fe/H]. First, higher metallicity spectra are more sensi- 
tive to errors in T^s and log g (see § 16. Sp . Second, higher 
metallicity spectra exhibit more complex absorption from 
molecular transitions not seen at lower metallicity. 

In order to quantify the total error on [Fc/Hjjyjpg, for 
each cluster we find the systematic error Jsys[Fe/H]j^j^g 
that satisfies 
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Fig. 9. — ["^/Fe] ]y[ jj^3 vs. individual a element abundances determined from HRS. The dotted lines are one-to-one. The symbols, colors, 
and references are the same as in Fig. [3] 



[Fc/H]mrs - ([Fc/H] 



MRS/ 



'((JfitiFe/H] 



MRS 



f + (<5,ys[Fe/H]^j,s)2 



(10) 

Figure [13] shows the results. (5sys[Fe/H]j^j^g appears to 
approach 0.1 at low [Fe/H] and increase as [Fe/H] rises 
above —1.6. We choose to fit a function that asymptoti- 
cally approaches a constant at low [Fe/H] and increases 
linearly at large [Fe/H]. A function that satisfies these 



requirements is In (e + i 
with uniform weighting is 



Vbx\ 



<5sy.[Fe/H]Mj^g -- 
The total error is 



In e 



,0.1 



The least-squares fit 
(11) 



-0.758+1.86 [Fc/H] 



'5tot[Fe/H] 



MRS 



(<5fit[Fe/H]^j,s)V(4ys[Fe/H]- 



JMRs) • 

(12) 
The error bars on [Fe/H]j^j^g for all plots and tables in 
this article are calculated from Eas.fTT]and[T2l 



6.2. Total Error on [a/Fe] 

We repeat the same procedure to find the systematic 
error in [c«/Fe]]yjj^g. However, we compute the deviation 
from the [Q!/Fe]jjf,g — accounting for HRS error — rather 
than the mean of [«/Fe]]yjj^g. 



[a/Fe] 



MRS 



[a/Fe] 



HRS 



'((5fit[a/Fe]]^j^g)^ + (4ys[a/Fe]j^j^s)^ + ((5[a/Fc]jjRs)^ 



(13) 



The value that satisfies this equation is 5sys[a/Fe]j^j^g 
0.03. The total error is 



MRS 



((5fit[a/Fe]jy„g) -f (<5sys[a/Fe]]^Rg) 



(5tot [a/Fe] 

(14) 
The error bars on [Q^/Fe]j^p{g for all plots and tables in 
this article are calculated from Eq. [TH 

6.3. Errors from Atmospheric Parameters 

Errors in Toff and log g have the potential to change the 
measured abundances significantly. The total [Fe/H]j^j^g 
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[Mg/Fe], [Si/Fe], [Ca/FeL^and [Ti/Fe] HRS measurements for indi- 
vidual stars. Eas.l8landl9lgive the formulas for deriving [a/Fe]jjjjg 
and its error from the individual element measurements. The dot- 
ted line is one-to-one. T he s ymbols, colors, and references are the 
same as in Fig. |3] See S IS. 31 for a discussion of the large values of 
[a/Fe]jjj^g for M15 (green points). 
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Fig. 12.— Mean cluster values of [a/Fe] vs. [Fe/H] for MRS 
(black) and HRS (pink). For comparison, we show the relation for 
MW stars from Venn et al. (2004). The dashed line is the average 
[a/Fe] in a moving window with a half width of 0.5 dex in [Fe/H], 
and the shaded region is the rms spread. [cf/Fe]]y[j5g is very near 
the MW halo plateau value of [a/Fe] = 4-0.3 for every cluster. 




Fig. 11. 
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The range of [a/Fe] in these GCs is very small, but the points are 
roughly consistent with the one-to-one line (dotted line). See §[33] 
for a discussion of the large value of ( [a/Fe] jjp^g ) for M15. 

error estimates in iJ IG.ll account for random error about 
the true values of Tcff and logg, but not systematic off- 
sets. The total [a/Fe]jyjj^g error estimates in ij |6.2l do 
account for both random and systematic error in the at- 
mospheric parameters to the extend that Teff and logg 
in the comparison HRS studies are accurate. 

We quantify the abundance error introduced by errors 
in atmospheric parameters by varying TcH' and logg for 
all stars in the GC sample. We recompute abundances at 
Teff ±125 K, Teff ±250 K, log 5 ±0.3, and log g± 0.6. Ta- 
ble [6] shows the differences in [Fe/H] and [a/Fe] between 
the altered and unaltered atmospheres. The numbers 
presented are the mean difference and standard devia- 
tions for all stars in all seven GCs. 

As expected, increasing Teff causes [Fe/H] to increase 
because the synthetic atmosphere must have a higher 
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Fig. 13. — The systematic error in [Fe/H] as a function of the 
mean [Fe/H] for each cluster. The line is the best fit of the function 



TABLE 6 
Errors from Atmospheric Par.ameters 



Atmospheric Error 



<5[Fe/H] 



5 [a/Fe] 



Teff - 250 K 
Teff - 125 K 
T,fi + 125 K 
Teff + 250 K 
logg — 0.6 
logg-0.3 
logs + 0.3 
logg + 0.6 



-0.21 ±0.08 
-0.10 ±0.04 
+0.11 ±0.05 
+0.20 + 0.08 
-0.06 ± 0.06 
-0.03 ± 0.04 
+0.04 + 0.03 
+0.07 + 0.07 



+0.01 + 0.12 
+0.01 + 0.08 
-0.03 + 0.06 
-0.07 + 0.10 
+0.00 + 0.07 
+0.01 + 0.06 
-0.01 + 0.05 
-0.01 + 0.07 



density of absorbers to compensate for the temperature- 
induced weakening in line strength. Increasing logg also 
causes [Fe/H] to increase because a higher electron pres- 
sure increases the density of H~ ions. H^ is the dominant 
source of continuous optical and near-infrared opacity in 
these cool giants. Therefore, as logg increases, the de- 
creasing ratio of line opacity to continuous opacity de- 
presses line strength. The synthetic atmosphere needs 
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Fig. 14. — Examples of the Ca II triplet region of spectra 
at three different signal-to-noise ratios. These are three different 
stars in NGC 2419, obtained with 20 minutes of exposure time on 
DEIMOS. Shown in each panel is the S/N, / magnitude, and mea- 
sured [Fe/H]. The value fro m one star (not show n) observed with 
HRS is [Fc/H]jjRs = -2-32 (IShetrone etalll200ll ') . 

to be more metal-rich to compensate. We also point out 
a few trends with Toff, logg, and [Fe/H]: 

1. (5[Fe/H]j^j^g is gradually less sensitive to ST^g as 
Toff increases. 

2. (5[Fe/H]j^pj^g does not show a trend with S\ogg as 
logg increases. 

3. (5[Fe/H]j^pjg does not show a trend with ST^s as 
[Fe/H] increases. 

4. (5[Fe/H],^j^g is gradually more sensitive to Slogg as 
[Fe/H] increases. 

5. On average, ^[Q!/Fe]]yjj^g has the same sign as ST^g 
for Toff < 4600 K and the opposite sign otherwise. 
(5[a/Fe]]y[j^g is most sensitive to STcs at low Tcff. 

6. '5[a/Fe]]yjj^g does not show a trend with Slogg as 
logg increases. 

7. (5[a/Fe]]y[pjg does not show any trends with [Fe/H]. 

6.4. Effect of Noise 

Stars at the tip of the RGB in MW GCs are easy tar- 
gets for high resolution spectrometers. The intended tar- 
gets of the MRS method are much fainter. Therefore, we 
explore the effect of noise on the measurement of [Fe/H] 
and [a/Fe]. To estimate S/N, we compute the absolute 
deviation from 1.0 of all pixels in the continuum regions 
(see ij |3.4p . We clip pixels that exceed 3 times this mean 
deviation. The S/N per pixel is the inverse of the mean 
absolute deviation of the remaining pixels. To convert to 
S/N per A, multiply by 1.74, the inverse square root of 
the pixel scale. 

Figure [14] shows the Ca II triplet region of three 
DEIMOS spectra of stars in NGC 2419 in three differ- 
ent S/N regimes. The three strongest lines are the Ca II 
triplet, and the five weaker lines are individual or blended 
Fe I transitions. We show the Ca II triplet to demon- 
strate spectral quality, but we do not use any spectral 
information from any of the three lines because we do 
not model them accurately. Instead, we use the weaker 
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Fig. 15. — The difference between individual measurements of 
[Fe/H]j^jj^g and the cluster's mean ( [Fe/H] j^^ j^g ) vs. signal-to-noise 
ratio. The dashed line is the rms in a moving window with a half 
width of 50 A . The dotted line is the average error determined 
from Eas. llll and ll2l in the same moving window. 
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Fig. 16. — Same as Fig. 1151 for measurements injected with 
artificial noise. The ranges of both axes are different from Fig. 1151 
The dashed line is the rms in a moving window with a half width 
of 5 A . The dotted line is the average error determined from 
Eas. llll and ll2l in the same moving window. 



metal lines, any one of which is difficult to identify in low 
S/N spectra. As an ensemble over ^ 2800 A of spectral 
range, they provide accurate abundances. Below, we dis- 



in which even the Ca II 



cuss spectra at S/N ^^ 10 A 
triplet is barely identifiable. 

Figure [T51 shows the difference between [Fe/H]j^j^g for 
an individual star and the mean [Fe/H]jy[j^g that we mea- 
sure for its cluster (A[Fe/H]]y[pg) vs. S/N. We also in- 
clude the rms of all points in a moving window with 

a hah width of 50 A" \ At S/N > 200 A^\ the rms 
scatter in A[Fe/H] is about 0.15 dex. The rms rises to 

about 0.20 dex in the range 20 A" ^ < S/N < 200 A"\ 
Errors such as these are about the magnitude of typi- 
cal HRS abundance errors. We also include the average 
<5tot[Fe/H]]y[j^g in the same moving window. The close 
agreement between the two error averages demonstrates 
that we determine (5tot[Fe/H]j^jjg well. 
In order to better scrutinize noisier spectra at S/N < 

20 A , we inject the raw spectra with Gaussian ran- 
dom noise proportional to the square root of the mea- 
sured variance in each pixel. On average, the S/N of a 
given spectrum decreases by a factor of 10. We repeat 
the complete analysis, including continuum determina- 
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Fig. 17. — Differences between [o/FeJi^jjj^g and [o/Fe]jjjj^g vs. 
signal-to-noise ratio. The error bars are the sum in quadrature of 
i5tot[«/Fe],^jj^g ("Eq. I14I I and 5[«/Fe]jjjjg (Eq.[9ll- The dashed line is 
the average of -^ ( [a/Fe] j^j j^g — [a/Fe]] 
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moving window with a half width of 0.5 dex in S/N. At S/N > 320, 
5[a/Fe]jjjj^g completely accounts for the total error, and the argu- 
ment of the square root becomes negative. The symbols, colors, 
and references are the same as in Fig.[3]except that faded symbols 
represent measurements from spectra injected with artificial noise. 

tion and velocity cross-correlation. Figure [TBI shows that 
rms scatter gradually increases as S/N decreases, but 

even at S/N ~ 10 A , the rms scatter is only 0.5 dex. 
The Ca II triplet is barely identifiable, and the analy- 
sis includes only lines weaker than the triplet, yet the 
spectra still contain enough information to yield decent 
metallicity estimates. 

We also investigate the effect of noise in estimating 
[a/Fe]. We repeat the comparison in Fig. [TO] in which 
the sample is restricted to those stars also observed 
with HRS. In Fig. [TTj we plot the difference between 
[a/Fe]]y[j^g and [a/Fc]jjj^g vs. S/N. In the same figure, 
we also plot the results from the same stars injected with 
artificial noise (faded symbols). The dashed line is the 
scatter about zero in a moving window with a half width 
of 0.5 dex in S/N, correcting for the error on [a/FeJjjpg: 

\/([a/Fe]MRS - ["/FelHRs)' " ('5[a/Fe]HRs)'- At 

S/N > 320, (^[a/FeJjjpg becomes comparable to 

-1 



I [a/Fe] 



MRS . 



a 



/Fc] 



HRS I 



Even at S/N 



10 



(the typical S/N of a one-hour exposure of an RGB 
star in M31), the average error is ~ 0.2 dex, small 
enough to distinguish between the halo plateau value of 
[a/Fe] ~ +0.3 and [a/Fe] - 0.0. 

A Gaussian noise model is a simplified representation 
of the actual noise model. In reality, the spectral er- 
ror contains a systematic component, particularly around 
night sky lines. A thorough noise model is beyond the 
scope of this article, but we intend to address this issue in 
future work. In the meantime, the Gaussian noise model 
gives a good approximation of the abundance precision 
we expect with lower S/N spectra. 

7. APPLICATIONS 

GCs are already well studied at high resolution. MRS 
abundance measurement is not intended to provide more 
accurate abundances for GCs. After all, it has only two 
dimensions of abundance, although future versions may 
even allow estimates of individual element enhancements, 
such as Al, O, Na, Mg, Ca, Ti, La, Ba, and Eu. However, 
the real strength of MRS is to probe large samples at 



large distances. 

Upcoming studies will focus on the metallicities and 
a enhancements of the dwarf galaxies of the MW and 
M31. We also intend to explore the differences between 
the chemical properties of different kinematic compo- 
nents of M31: the cold dSphs, the disrupted satellites 
and streams, and the hot halo. We intend for the MRS 
measurements to be a direc t test of che mical evolution 
simulations, such as those of I Font et al.l (|2006l). 

In exploring these diverse systems, we will need to rec- 
ognize some shortcomings of the MRS technique. First, 
we do not model TiO absorption, and we have discarded 
all stars that show TiO. In future applications, we intend 
to continue discarding stars with Toff < 4000 K due to 
the complexity of modeling TiO absorption. Second, the 
errors at the high er metallicity of M71 ([Fe/H] — —0.76, 
iPritzl et"aLl I2005D become somewhat large (0.2 dex), 
though not larger than typical errors for photometric or 
Ca II triplet-based metallicities. Both problems will af- 
fect metal-rich sample s such as the in ner halo of M31 
(([Fe/H]) = -0.47. :Kalirai et all [2001 . Finally, sam- 
ples of very faint stars may require spectral coaddition 
in bins oiTes and log 5 to achieve S/N ratios high enough 
to make measurements of a enhancement with precision 
greater than 0.3 dex. Although enhancements of individ- 
ual stars would be unrecoverable, coadded spectra would 
contain enough information to estimate a stellar popula- 
tion's overall chemical properties. 

We have demonstrated that metallicity measurements 
of medium resolution spectra via spectral modeling ap- 
proach the accuracy and precision of measurements from 
HRS. Furthermore, the leverage of a large number of ab- 
sorption lines from many different elements allows abun- 
dance measurement in multiple dimensions. Presently, 
the method estimates a enhancement, but future ver- 
sions may provide individual element enhancements. 
Spectral modeling is superior to spectrophotometric in- 
dices, such as the Ca II triplet EW. First, the preci- 
sion from modeling is comparable or better at the same 
S/N. Second, modeling requires no empirical calibration. 
Therefore, the result is not tied to intrinsic parameters 
(e.g., [Ca/Fe]) of the calibrators. Additionally, modeling 
is not subject to the range of metallicities of the cali- 
brators. Whereas empirical calibrations to GCs will be 
inaccurate at [Fe/H] < —2.4, the MRS method is appli- 
cable to arbitrarily low metallicities, whi ch may exist in 
some of the lowest mass MW satellites (jSimon fc Gehal 
|2007| ). These advantages make medium resolution spec- 
tral modeling the only abundance technique that will be 
able to test theories of hierarchical structure formation 
through large samples of individual stars at large dis- 
tances. 
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TABLE 7 
Abundances for Individual Stars 



a (2000.0) (5 (2000.0) 



_ga 



[Fe/H], 



[a/Fe]- 



[Fe/H]j 



[a/Fe], 



HRS reference 



M13 



16''40™52.1= 


+36°29'25" 


17.31 ± 0.04 


16.64 + 0.03 


15.82 + 0.01 


5336 


3.09 


-1.69 + 0.13 


+0.33 + 0.07 








16'^40""58.4= 


+36°26'04" 


16.70 + 0.02 


16.03 + 0.01 


15.17 + 0.01 


5204 


2.80 


-1.86 + 0.11 


+0.41 + 0.05 








16'^41"'01.6= 


+36°28'00" 


17.37 + 0.04 


16.79 + 0.02 


15.95 + 0.01 


5260 


3.13 


-1.71 + 0.12 


+0.28 + 0.07 








16'^41'^02.6= 


+36° 26' 16" 


15.68 + 0.01 


14.92 + 0.01 




4965 


2.26 


-1.65 + 0.12 


+0.22 + 0.04 


-1.44 + 0.05 


+0.22 + 0.03 


Cohen & Melendez 


16li4imQ4gs 


+36°27'45" 


17.48 + 0.04 


16.92 + 0.01 


16.10 + 0.01 


5318 


3.20 


-1.82 + 0.12 


+0.23 + 0.09 








16'i41'"05.8= 


+36°26'30" 


17.47 + 0.04 


16.86 + 0.01 


16.05 + 0.01 


5379 


3.19 


-1.68 + 0.13 


+0.27 + 0.07 








16'^41™06.2= 


+36°25'22" 


15.30 + 0.01 


14.46 + 0.01 




4829 


2.02 


-1.65 + 0.12 


+0.23 + 0.04 


-1.44 + 0.05 


+0.17 + 0.04 


Cohen & Melendez 


16'^41""06.5= 


+36° 28' 14" 


14.29 + 0.01 


13.25 + 0.01 




4522 


1.37 


-1.65 + 0.12 


+0.11 + 0.04 


-1.63 + 0.06 




Sneden et al. (2004' 


16h4imo9.7s 


+36°26'45" 


15.52 + 0.02 


14.73 + 0.01 


13.79 + 0.01 


4997 


2.19 


-1.67 + 0.12 


+0.30 + 0.04 








161141^09.9= 


+36°27'42" 


17.31 + 0.03 


16.73 + 0.02 


15.91 + 0.01 


5352 


3.13 


-1.52 + 0.14 


+0.13 + 0.08 









Note. — Tablc[7]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding form and content. 
^ P. B. Stetson has generously provided this photometry, which is preliminary pending his own publication of these data. See ^ |5.3| for a discussion on [a/Fe]pjDg measurements of lSneden et aU 
(fT997h . fSneden et al.l ||2000D . and lSneden et al.l (f200l^ . 






